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WALL-PRESSURE-ARRAY MEASUREMENTS BENEATH 
A SEP ARATNG/REATT ACHING FLOW REGION 


Laura M. Hudy*, Ahmed M. Naguib^ 

Michigan State University 
East Lansing, MI 48823 

William M. Humphreys Jr. { , Scott M. Bartrarh 1 
NASA Langley Research Center 
Hampton, VA 23681 

A database of wall-pressure array measurements was compiled for studying the space- 
time character of the surface-pressure field within a separating/reattaching flow region. The 
experimental setup consisted of a long splitter plate instrumented with an array of 80 flush- 
mounted microphones located within the wake of a fence. Data were acquired for a Reynolds 
number of 7885, based on the fence height. Two distinctive regions, defined based on their 
location relative to the position of the mean reattachment point (x r ) of the shear layer, emerged 
from this investigation. Upstream, from the fence to !4x r , the surface-pressure signature was 
dominated by large time scale disturbances and an upstream convecting velocity of 0.2 1U,,. 
Beyond !4x r , turbulent structures with small time scales and a downstream convection velocity of 
0.57U„ generated most of the pressure fluctuations. There was evidence that these structures 
began to form around V$x r and grew in strength and size with downstream distance before 
reattaching on the plate. Only the time-averaged results from the microphones have been 
examined hitherto and will be presented. 


Introduction 

Separating/reattaching flows produce large 
pressure fluctuations on the underlying surface. These 
fluctuations could cause significant vibration of the 
surface and subsequent generation of noise. To predict 
and/or control such vibration and noise effects one 
needs to understand the spatio-temporal character of the 
surface-pressure field. The present study addresses this 
issue through the use of a wall-microphone array to 
resolve the surface-pressure field both spatially and 
temporally in a basic separating/reattaching flow 
geometry. This provides further contribution to the 
bulk of the literature in this area, which has been 
primarily based on only one- or two-point 
measurements. 

The flow geometry investigated in the current 
study consists of a splitter-plate attached to, and 
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downstream of, a fence that is perpendicular to the 
flow, as shown in Figure 1. This model was used 
because it has a separation bubble that is elongated in 
the streamwise (x) direction, and hence the 
development of the wall-pressure field within the 
bubble can be resolved properly using a sensor array 
with inter-sensor spacing that is not too small to realize. 
Also in this geometry the boundary layer thickness is 
much smaller than the fence ('step') height, and 
therefore the boundary layer details have minimal affect 
on the flow field. 



Figure 1. Splitter plate with fence flow geometry 
investigated and ideal two-dimensional flow field 

Cherry et al. 1 made two-point unsteady surface 
pressure measurements in a separating/ reattaching flow 
region. Their test model geometry was a blunt-face 
splitter plate, which is similar in nature to the splitter- 
plate-with-fence geometry. The Reynolds number used 
in the experiment was based on the thickness of the 
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splitter plate and was held within the range of 3.2x1 0 4 ± 
0.2xl0 2 . Cherry et al } demonstrated the convective 
nature of the surface-pressure imprint associated with 
the downstream motion of the shear layer structures 
through cross-correlation analysis between two signals 
from microphones spaced apart in the streamwise 
direction. One microphone was fixed near the 
reattachment and the other was moved to five different 
positions from upstream to downstream of the fixed 
microphone. The peak in the cross-correlation function 
was seen to shift to a smaller time delay as the movable 
microphone was traversed in the downstream direction, 
indicating the convective nature of the flow structure 
dominating the wall-pressure generation process. 
Cherry et al? also observed low- and high-frequency 
peaks in the power spectrum measured close to 
separation and near reattachment, respectively. They 
attributed the low-frequency signature to the flapping of 
the shear layer associated with the growth and decay of 
the separation bubble. 

Farabee and Casarella' studied the fluctuating 
wall pressures in a forward- and backward-facing step 
flow using a flush-mounted B&K 1/8 in. condenser 
microphone. Using frequency domain analysis, 
Farabee and Casarelkr described the characteristics of 
the wall-pressure field as variable with x distance along 
the wall. Close to separation, the spectra showed the 
highest level of energy at low frequencies; whereas, 
farther downstream the spectrum containing the largest 
energy was found at reattachment. This was a 

manifestation of the increase in the energy of the 
organized, turbulent structures as the flow convected 
downstream. A corresponding shift was seen in the 
spectrum as the dominance of the low-frequency 
disturbances gave way to the dominance of the high- 
frequency structures downstream. A convection 
analysis at different positions along the model 
downstream of the step showed that the pressure 
fluctuations close to separation were associated with the 
re-circulating low-speed fluid and not the high-speed 
fluid in the shear layer. However, Farabee and 
Casarella" commented that the convection velocity was 
always in the downstream direction, indicating that the 
pressure fluctuations were not originating from the 
reverse flow within the re-circulating bubble. 

Driver et al. 3 , in their backward-facing step 
study of the time-dependent character of the separated 
shear layer at a Reynolds number of 37,000, based on 
the step height, noticed abnormal contraction and 
elongation of the separation bubble due to the 
shortening and lengthening of the reattachment length. 
This was labeled as the flapping motion of the shear 
layer with amplitude estimated to be 20% of the shear 


layer width. They used surface pressure measurements 
along with velocity measurements to show that there 
was a definite low frequency disturbance associated 
with the shear layer flapping, but that it contributed 
very little energy to the overall pressure fluctuations. 

Heenan and Morrison 4 investigated wall- 
pressure fluctuations behind a rearward-facing step and 
passive control of these fluctuations using a permeable 
surface (at Reynolds number equal to 1.9xl0 5 based on 
the step height). Heenan and Morrison 4 found an 
upstream convection velocity close to separation using 
phase-angle analysis. They identified negative phase 
angles (with respect to a microphone signal measured 
immediately behind the step) at low frequencies and at 
locations from separation up to 0.4x r in the 
impermeable case. This is the only study found to date 
that describes an upstream convecting velocity. 

Lee and Sung 5 used a 3 2 -microphone array 
downstream of a backward-facing step to measure wall- 
pressure fluctuations in the streamwise and spanwise 
directions. Spatio-temporal statistics were completed 
on this comprehensive data set for a Reynolds number 
of 33,000, based on the step height. Lee and Sung 5 
observed the same phenomenon in their experiment as 
experienced by earlier investigators of backward- 
facing-step studies. The RMS pressure fluctuations 
rose sharply starting around 0.5x r and peaked in the 
vicinity of reattachment, decaying beyond that point. 
Pressure spectra revealed low-frequency dominance 
close to separation, presumably due to the flapping of 
the shear layer. Farther downstream, the spectra were 
dominated by high-frequency components. In terms of 
convection velocity, Lee and Sung 5 calculated a 
downstream convection velocity of 0.6U^ at high 
frequencies and they found no evidence of an upstream 
convection velocity. Although from their phase plot, 
(used to determine the convection velocity) there were 
many singularities (phase discontinuities) at low 
frequencies. This was not the case at the higher 
frequencies. 

The characteristics described thus far in these 
surface pressure measurement studies have also been 
identified by authors that used different measurement 
techniques within similar type of geometries. In these 
separating/reattaching flow studies the measuring 
techniques used included hot-wire and pulsed-wire 
anemometry, skin-friction measurements, and particle 
image velocimetry. Castro and Haque 6 , Eaton and 
Johnston , and Spazzini et a I s are a few of the authors 
who also observed very large-scale, low-frequency 
motion close to separation and small-scale, high 
frequency motion close to reattachment using these 
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techniques. In all studies, the low-frequency signature 
was attributed to the flapping of the shear layer. 


Experimental Set-up 

The present experiment was completed in the 
open-circuit. Subsonic Basic Research (Wind) Tunnel 
(SBRT) at NASA Langley Research Center in 
Hampton, Virginia. This low-speed wind tunnel has a 
6:1 contraction section that measures 3.54 m in length. 
Located downstream of the contraction is a 0.57 m wide 
by 0.82 m high by 1.85 m long test section. Flow 
entering the wind tunnel is straightened using 
aluminum honeycomb flow straighteners and a double 
row of wire mesh turbulence screens. The fan is driven 
by a 200hp motor enabling flow speeds in the wind 
tunnel up to 60m/s. For this study, the flow speed (U„) 
used was 15m/s, resulting in a Reynolds number of 
7885, based on the step height of the model. The 
corresponding turbulence intensity is less than 3%. 

A schematic of the test model is shown in 
Figure 2, where x represents the streamwise distance 
measured from the fence, y is the normal distance from 
the splitter plate, z is the spanwise distance, h f 
represents the step height above the splitter plate (8 
mm), and 2H represents the total fence height (35 mm). 


Flow 

— ► 



X 


1 


Figure 2. Splitter-plate-with-fence configuration 


The design of the model was symmetric with 
respect to top and bottom. It was constructed out of 
aluminum using a 12.7 mm thick aluminum skeleton 
support covered with 3.175 mm thick aluminum sheets, 
which provided storage space for the microphone 
wiring and static pressure tubing inside the model. The 
total length of the model was 160h f or 73H and its 
width was 44h f or 10H. Endplates were placed on the 
sides of the splitter plate to improve the two- 
dimensionality of the flow, according to Castro and 
Haque 6 , resulting in a model aspect ratio of 36. The 
blockage ratio, as defined by Smits 9 who found that the 
reattachment distance decreased with increasing 
blockage ratio, was around 2%. This resulted in a 
reattachment length of approximately 0.2 m (25.6hf), 
which ensured that the reattachment point would be 
contained within the extent of the microphone array. 

The splitter plate was outfitted with 80 
microphones and 80 static pressure taps. The 
configuration and numbering of the microphones and 
the static pressure taps are shown in Figure 3. 


The microphone array consisted of Panasonic 
(WM-60A) omnidirectional back electret condenser 
microphone cartridges with a nominal sensitivity of - 
44±5 dB (relative to lV/Pa) and a bandwidth of 20- 
20,000 Hz. The microphones, each with a sensing 
diameter of 2mm, were flush mounted and were used to 
record the fluctuating pressure on the surface of the 
plate. The center row consisted of 28 microphones 
spaced 1.2h f apart center to center in the streamwise 
direction, starting at 0.6h f downstream of the fence. On 
either side of the centerline, there were 2 rows each 
containing 13 microphones. These 13 microphones 
were spaced 2.4h f apart center to center, except for 
three microphones in each row that were spaced 1.2hf 
apart center to center. The spanwise spacing between 
the five rows was 6.4h f . 


Upstream (fence location) 
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Figure 3. Top view of the instrument plate layout 
consisting of microphones and static pressure taps 

Each microphone was calibrated against a Vi" 
B&K microphone with known sensitivity before being 
placed in the model. The outputs of sixty-four of the 
microphones were connected to eight National 
Instruments A/D Boards (SCXI 1141), placed in a 
SCXI 1001 chassis. Each board had an input signal 
range of ±5V, eight differential analog-input channels, 
and a variable channel gain that was set to one for this 
experiment. The highest sampling rate the board was 
capable of was 1.25 MHz for one channel. For the 
sixty-four channels, the corresponding maximum 
sampling rate was 19,53 1 Hz per channel resulting in an 
average time delay between successive channels of 2.6 
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|is. This delay, however, was small compared to a 
characteristic time scale of the flow eonveetion over a 
distance equal to the microphone spacing: A/LL = 635 
|is; where A is the spacing between successive 
microphones in the streamwise direction. Data for the 
experiment was sampled at 12kHz for 15 seconds, with 
the cut-off frequency anti-aliasing filters set at 5kHz. 

The splitter plate was also instrumented with 
40 static pressure taps on the top and bottom for a total 
of 80 taps. Static pressure measurements were 
primarily used to align the model in the tunnel and 
estimate the reattachment length. The location of the 
taps, which mirrored each other on top and bottom, is 
depicted in Figure 3. The static pressure taps were 
coupled to a 48-port Scanivalve that was connected to a 
Setra 239 series pressure transducer. The transducer 
measured differential pressure in the range of 0- 
25.4mm H 2 Q, outputting a corresponding 0-5V signal. 


Although these studies were conducted at 
different Reynolds numbers and had various blockage 
ratios, resulting in a difference in the reattachment 
lengths, the pressure distribution for all experiments 
correlate well when plotted using the C p * coefficient, 
even in the case of the blunt-face-splitter -plate 
geometry. Thus, a universal pressure coefficient (C p *) 
value may be found at the mean reattachment location 
(x/x r =1), within the band of data scatter (about ± 5%). 
This value was determined to be approximately C p * = 
0.35. By applying this C p * value to the present data, 
the reattachment distance could be determined within ± 
5% uncertainty as shown in Figure 4. Therefore, the 
reattachment distance was determined to be roughly 
205mm or 25.6h f . Static pressure tap and microphone 
#21 is the port closest to this reattachment value and 
thus was used in the normalization of the present data 
throughout the study. 


Results and Discussion 

Reattaclunent length 

The reattachment length (x r ) is an important 
parameter for the present flow geometry and, as shown 
by Ruderich and Femholz 10 , appears to be the 
appropriate length scale for this flow field. Therefore, 
before analyzing the data it was necessary to estimate 
the reattachment length for the purpose of normalizing 
the results. This was done utilizing the pressure 

coefficient (Cp*) used by Ruderich and Femholz 1 " in 
their presentation of C p results. The form for collapsing 
of C p * was first proposed by Roshko and Lau". As 
Ruderich and Fernholz 10 explain, the mean pressure 
distribution results from different long-separation- 
bubble studies: 

* C p- C p. mi n 

v i - r 

1 ^p, min 

where C p is defined as (prPr)A%pLL 2 ), p s is the surface 
pressure along the model, p r is a reference pressure, 
measured with a static pressure tap located at the exit of 
the contraction, and C pjtlin is the minimum C p in the 
mean pressure distribution. Figure 4 shows the C p * 
distribution for seven different studies, including the 
present one, as a function of the distance along the 
splitter plate normalized by the reattachment length. 
All seven studies, except Cherry et al. 1 , used a splitter- 
plate-with-fence configuration. Cherry et al } 
investigated a blunt-face splitter plate, with the 
thickness of the plate denoted by D. The legend 
displays details about each of the studies for which data 
are shown in Figure 4. 
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Figure 4. C p * distribution compared with six different 
studies 


Mean and Fluctuating Pressure Distributions 

Figure 5 shows the mean pressure distribution 
for three different studies, including the present 
measurements and data from Castro and Haque 6 and 
Cherry et al ., as a function of the distance along the 
instrument plate normalized by the reattachment 
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distance. There is good agreement between the current 
results and the Castro and Haque 6 data, which indicates 
the consistency of the mean-flow behavior around the 
constructed model with earlier studies of the same 
geometry. There is some offset between the present 
data and that from Cherry et al . 1 This offset could be 
due to the difference in model geometry. Comparison 
with Cherry et al. 1 , however, was important because it 
was the only detailed study found with more than one- 
point unsteady surface pressure measurements in a 
separating/reattaching flow geometry similar in nature 
to the present experiment (i.e., one where the boundary 
layer thickness at separation is much smaller than the 
step height). Cherry et al. 1 conducted two-point 
measurements at different spacings on a splitter plate. 
Kiya and Sasaki 1 ' also studied the flow state over a 
blunt-face splitter plate using extensive single- and two- 
point measurements of surface-pressure and velocity 
measurements. However, they displayed most of their 
data in velocity and velocity-pressure correlation plots. 
It should also be noted that although the study of Lee 
and Sung 3 , discussed in the Introduction, did utilize a 
32-microphone array, their flow geometry was a 
backward facing step and not a splitter plate, with or 
without a fence. 
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Figure 5. Streamwise distribution of the mean pressure 
coefficient for the current study compared to Castro and 
Haque 6 and Cherry et al. 1 

The root-mean-square (RMS) of the pressure 
fluctuation is shown in Figure 6 for the present study 
and that of Cherry et al . 1 The RMS data are plotted as 
Cpi— P'mu/l^pU - in Figure 6. The present data exhibit 
some scatter around the trend line, which is believed to 
be associated with the uncertainty of the microphone 
calibration procedure used to determine the sensitivities 
of the microphones. This uncertainty was found to be 
approximately 7%. There is qualitatively good 
agreement between the two data sets, although there is 
some difference as seen in the static pressure 


measurements. This offset is largest in the vicinity of 
the peak in the C p > values and could be due to the 
difference in the model geometries selected for each 
study. 

At the point of separation, the shear layer is 
laminar and relatively far away from the splitter plate 
wall-pressure sensors. At this location, the 
microphones detect low RMS pressure fluctuations 
directly behind the fence. It is unknown what causes 
these pressure fluctuations, but it has been theorized in 
the literature that the unsteadiness, or ‘flapping’, of the 
shear layer may in fact produce some of the wall- 
pressure activity seen in this region. The region 
referred to here is the distance from the fence up to 
about a quarter of the reattachment distance (!4x r ), 
where the RMS values are relatively flat for both data 
sets shown in the graph. This region is significant and 
will be referred to in the upcoming sections. 
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Figure 6. Streamwise distribution of the coefficient of 
RMS pressure fluctuations 

Beyond the l Ax T , there is a rise in the RMS 
pressure fluctuations, as seen in both data sets in Figure 
6. It is believed that in this region the surface pressure 
fluctuations are predominately associated with the 
vortical structures of the separated shear-layer. As 
these structures convect downstream, growing in size 
and strength and moving closer to the wall, they 
produce an increasingly strong wall-pressure signature. 
The signature reaches a maximum level in the vicinity 
of where the flow impinges, or reattaches, on the plate 
as described by Farabee and Casarella - . It is well 
documented that the peak RMS value occurs slightly 
upstream of the reattachment point in both splitter -plate 
and backward-facing-step studies. Heenan and 
Morrison 4 , in their backward-facing-step study, found 
the maximum RMS value to occur approximately one- 
step height upstream of reattachment. Beyond 
reattachment, the RMS values decrease slowly as the 
energized structures from the shear layer decay and 
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diffuse downstream. The flow then takes on boundary 
layer characteristics once the shear layer reattaches. 

Auto -correlation and Integral Time Scales 

To identify the dominant time scales in the 
measured wall-pressure time records, an auto- 
correlation analysis was conducted. The gray-scale 
contour map in Figure 7 shows the auto-correlation 
coefficient (Rpy) for all 28 centerline microphones. 
The abscissa shows the distance along the splitter plate 
with respect to the reattachment point, while the 
ordinate shows the time shift normalized by the free 
stream velocity and the total fence height. The color 
bar indicates the values of the auto-correlation 
coefficient, which was obtained by normalizing the 
correlation function by the square of the RMS of the 
signal. The map makes it easy to see the variation of 
the width of the auto-correlation with downstream 
distance and the transition between long and short 
timescales. 



- 0.8 - 0.6 - 0.4 - 0.2 0 0.2 


(x-x )/x 
r r 

Figure 7. Contour map of the auto-correlation 
coefficient for all 28 microphones along the centerline 
of the model 

In the region immediately behind the fence, 
the auto-correlation function extent is wide and changes 
very little up to a distance of about 0.2 - 0.25x r (x-x,/x r 
= -0.80 to -0.75) behind the fence. Farther 
downstream, this width narrows significantly over a 
relatively short distance (roughly from 0.25x r to 0.7x r ; 
x-Xj/x,. = -0.75 to -0.30) as demonstrated by the 
focusing of the R p r p i contours towards XLL/2H = 0 line. 
Beyond this region, the contour lines remain 
approximately parallel to the constant x lines showing 
very little change in Rp. p . with additional increase in x. 
The region between x/x r = 0.25 to 0.7 roughly 
delineates the start and end locations of the change in 


the time scales of the flow structures dominating R p . p .. 
This region is referred to as the changeover region. 

The dominance of low-frequency disturbances 
directly behind the fence has been identified in a 
number of studies. These include Castro and Flaque 6 , 
Cherry et al. x . Driver et al. 3 , Eaton and Johnston 7 , 
Farabee and Casarella", and Lee and Sung 5 . Some of 
these studies have attributed these disturbances to the 
flapping of the shear layer as discussed in the 
introduction. Farther downstream, the organized shear 
layer structures grow in strength and move closer to the 
wall. These more energized structures impose a shorter 
time scale, than that encountered close to the fence, on 
the auto-correlation function. Thus, the increasing 
influence of these structures on the wall pressure 
appears to be responsible for the observed change in the 
Rp.p' within the changeover region. Past this region, in 
the vicinity of the reattachment location and farther 
downstream, the energy of the shear layer organized 
structures appears to saturate (as seen in the RMS plot 
in Figure 6). This is possibly why no substantial 
change in R p . p . is detected past the changeover zone. 

Figure 8 shows the streamwise variation of the 
integral time scale (x*) as a function of distance along 
the splitter plate with respect to the reattachment point. 
The integral time scale was derived from the auto- 
correlation function by finding the time at which the 
negative peak in the auto-correlation occurred with 
respect to the ordinate and multiplying that time value 
by two. Because of the even symmetry of the auto- 
correlation, x* gives the time delay between the two 
negative peaks in the auto-correlation function. An 
example of the method is shown in Figure 8a. The 
displayed error bars represent the uncertainty in 
locating the peak. A substantial error is encountered in 
the changeover region because the negative correlation 
peak in this region is flat, and therefore difficult to 
locate precisely in the presence of experimental data 
scatter (Figure 8b). 

Figure 8 is similar in nature to that produced 
by Castro and Flaque 6 using velocity measurements in 
the shear layer region. The results at the locations of 
the first four microphones show long time scales with 
relatively small error bars. This is the region closest to 
separation where other authors and this study have 
observed predominately low-frequency motion. Data 
from the next five microphones, associated with large 
error bars, are in the changeover area. The remaining 
microphones have shorter time scales with small error 
bars. This reaffirms the observations from the auto- 
correlation contour plot. That is, the signals from the 
microphones closest to separation are dominated by 
large time scales, whereas, the signals from the 
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microphones farther downstream are dominated by 
smaller time scales. In between, there is a changeover 
region starting roughly around l Ax r , which is the same 
region seen in the RMS statistics. This region is 
believed to be associated with the amplification and 
streamwise development of the vortical shear layer 
structures. 



Figure 8. Streamwise variation of the integral time 
scale 


Pow'er Density Spectra 

The power spectrum density at even' third 
microphone along the centerline is shown in Figure 9. 
This figure contains the power spectrum plots on a 
logarithmic scale for both axes. The ordinate is plotted 
relative to an arbitrary reference value and represents 
the power spectrum normalized by the square of iTpLL 2 
and the abscissa represents the frequency normalized by 
the fence height and the free stream velocity. The use 
of the arbitrary reference value for the ordinate provides 
a means by which many spectrum plots can be shown 
on the same graph without clutter. 

Figure 9 magnifies the spectra plots in order to 
get a closer look at the spectrum details and the shift in 
the peak frequency in the downstream direction. Close 
to separation, the peak in the frequency is seen to be 
roughly around f(2H)/LL = 0.02-0.03 (fx r /LL = 0.12- 
0.18) in the present study, which is the same peak seen 
by Cherry et al'. In the case of the backward-facing 
step, Lee and Sung 5 also found a similar peak frequency 
value close to separation at fH/LL = 0.015 (H = the 
height of the step), which when scaled by x r instead 
gives a value of fx r /LL = 0.1 1. This is comparable to 
the values given by Spazzini et al. 7 8 at fx r /LL, = 0.08 and 
Fleenan and Morrison 4 at fx r /LL = 0.1. Cherry et al . 1 
along with Fleenan and Morrison 4 , Driver et al. 3 , and 
Lee and Sung 5 , have associated this low-frequency peak 


with the flapping of the shear layer. Farabee and 
Casarella" suggested that the energy distribution in the 
spectra indicates that the wall pressure fluctuations 
close to separation were caused by the unsteadiness of 
the low-speed re-circulating flow, rather than the highly 
turbulent structures in the shear layer. This is 
consistent with the fact that these structures are only 
beginning to develop in this region and are most likely 
weak compared to the strength of the low frequency 
disturbance produced by the shear layer movement. 



covering the measurement range - streamwise direction 
is from top to bottom 

Farther downstream, the energy in the 
spectrum is located at higher frequencies as seen in the 
spectra in Figure 9. In particular, the peak in the 
spectrum occurs around f(2H)/U„ = 0.1-0.15 (fx/LL = 
0.6-0. 9), which is in agreement with the findings of 
Cherry et al. 1 , where the frequency is normalized using 
the total width of the splitter plate in their case. Lee 
and Sung 5 stated that the power spectrum reaches a 
maximum at fH/LL =0.068, or fx/LL = 0.5. Spazzini et 
a I s found their maximum at fx/LL = 1.0 along with 
Heenan and Morrison 4 . Driver et al. 3 recorded a peak 
value close to reattachment of fx r /LL = 0.6. This higher 
frequency peak has been attributed to the highly 
turbulent structures within the shear layer, as discussed 
previously. 

Furthermore, it is interesting to compare the x* 
values obtained from the auto-correlation to the spectral 
peak frequencies. This comparison is possible because 
of the inverse relationship between time and frequency. 
For large time scales, the x* value close to the fence 
was determined to be X*1L/(2I I) = 36.5; thus the 
corresponding frequency is f(2LI)/LL = 0.027, which is 
similar to the low-frequency peak found in the power 
spectra (f(2H)/LL = 0.02-0.03). Similarly, for the 
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smaller time scales farther downstream where the T* 
value is such that x*U 00 /(2H) = 7.3, resulting in a 
frequency of f(2I I)/LL = 0. 14 compared to the f(2H)/LL 
= 0. 1-0. 15 high-frequency peak from the power spectra. 
Castro and Haque 6 obtained velocity auto-correlation 
measurements at various positions in the shear layer 
using a pulsed-wire anemometer. From their results, 
they determined the time scale near separation to be 8 
when normalized with the reattachment distance. This 
value compares well with T*FL/x r = 6.2 calculated in 
the present study near the point of separation. The 
discrepancy could be due to the difference in the 
measuring techniques and/or the location where the 
data were recorded. 

Cross-correlation with respect to measurements at 
reattachment 

Figure 10 shows a gray-scale contour map 
(similar to Figure 7) of the cross-correlation coefficient 
of the signals from all 28 centerline microphones with 
that from the microphone closest to reattachment. 



- 0.8 - 0.6 - 0.4 - 0.2 0 0.2 


(x-x )/x 
r r 

Figure 10. Cross-correlation-coefficient color contour 
map for all 28 centerline microphones (reference 
microphone closest to reattachment) 

In the cross-correlation contour map, there is a 
main, positive correlation peak inclined at an angle and 
two negative lobes on either side of this main lobe. At 
each x location, the main peak is centered around time 
shift values corresponding to the largest positive 
correlation between the wall pressure signal at this x 
location and that at reattachment. The negative peaks 
give the time delay to the highest negative correlation. 
By finding the slope of the locus of the main positive 
(or negative) peak, an average downstream convective 
velocity can be calculated for the dominant turbulent 
structures regardless of their time scales (frequencies). 


The velocity is calculated from the slope of the locus 
line using 

1 m, (2FI) 

= — ^ — - ( 2 ) 

Uq UoqXj- 

where m s is the slope of the locus line determined from 
a least-squares line fit and U c is the convection velocity. 
Using equation (1), the convective velocity was 
determined to be 57% of the free stream velocity. 
Fleenan and Morrison 4 reported, using two flush- 
mounted pressure transducers in their backward-facing- 
step configuration, convective velocities ranging 
between 0.5U M and 0.6U«, depending on the position 
along the model. Lee and Sung 5 stated that convective 
velocities at high frequencies converged to a value of 
0.6U„ in their backward-facing-step study, using a 32- 
microphone array. Hwang et al. u estimated the 
average convective velocity in their blunt-faced flat 
plate to be approximately 0.5U„, based on flow 
visualizations using a high-speed camera with high 
framing rates. Cherry et al . 1 determined the convective 
velocity to be 0.5U^ from pressure-pressure cross- 
correlations. Generally, the convective values were 
cited to range from 0.5 to 0.6 of the free stream velocity 
in the literature, depending on model geometry, 
location of measurement, and measuring technique. 
However, consideration of the results cited above 
suggests that the convective velocity values of the 
dominant structures in the splitter plate/fence flow are 
more similar to the backward-facing step, than to the 
blunt-facing plate. 

Phase. Coherence, and more on Convection Speeds 

The cross-correlation results yield an average 
convection velocity associated with various time scales. 
In order to detennine the convection velocities 
associated with individual time scales (frequencies), the 
streamwise development of the phase angle (0) was 
examined for various frequencies. 0 at a given 
frequency and x location is computed from the cross 
spectrum between the microphone at that x location and 
the reference microphone (closest to the fence in the 
current results). However, before presenting the results 
for the phase plot, the range of frequencies for which 
reliable phase information can be calculated must be 
determined. This is accomplished by computing the 
coherence between the two signals for which the phase 
plot is to be obtained. The one-sided coherence 
between two signals obtained from N-point data records 
is defined as follows: 
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^p'i pS w — / 5 

^P'iP'i (k )^p' 2 p'2 Ck) (3) 

N 

k = 0 , 1 ,... — 

2 

where (|)p’ip '2 is the one-sided cross-spectrum between 
the two signals and <j)p’i p 'i and <j) P ' 2 P '2 are the 
corresponding one-sided power spectra. The coherence 
gives a measure of the ‘phase locking’ between two 
signals at a particular frequency. If the signals are 
perfectly correlated across all frequencies, then the 
coherence value will be unity over the entire range. In 
general, a coherence value larger than 50% at a certain 
frequency is indicative of the existence of a fairly well 
defined phase between the two signals at that 
frequency. Thus, for frequencies where r p 'i p -2 < 0.5, 
the calculated phase values may not be reliable and will 
generally be randomly scattered. Figure 1 1 shows the 
coherence for neighboring microphones (numbering 0 
and 1, 7 and 8, 16 and 17, and 26 and 27) to be high 
over a particular frequency range. 






Figure 1 1 . Coherence plot for four neighboring pairs of 
microphones 

The coherence plots in Figure 1 1 reveal that 
the signals measured by adjacent pairs of microphones 
are not coherent (phase-locked) across all frequencies. 
In fact, there is a sharp drop off in the coherence around 
f(2FI)/LL = 0.15, immediately behind the fence. This 
value increases with x up to f^FO/LL = 0.4 at x/x r = 
1.28. This drop-off is seen in all four coherence plots, 
which provides a good representation of the high- 
coherence frequency range found for any two 
neighboring microphones along the centerline. 
Therefore, the phase analysis will be conducted using 
neighboring microphones and will be constrained to the 


range up to f(2FI)/LL = 0.3 (this falls roughly in the 
middle of the range f(2Fl)/LL = 0.15-0.4). It should be 
noted here that this frequency range contains the low- 
and high-frequency peaks identified earlier in the power 
spectrum behind the fence and farther downstream, 
respectively. Thus, all the flow structures of interest to 
the surface-pressure generation process are contained 
within the frequency range where high coherence is 
seen for two neighboring microphones. 

Because of the use of pairs of neighboring 
microphones in calculating 0, there was no phase angle 
change greater than K between two signals at any 
frequency. This was confirmed by estimating a 
representative wavelength (X = U c /f) of the flow 
structures dominating the surface pressure 
measurements to be about 17hf, using the high- 
frequency peak f(2FI)/LL = 0.15 and the average 
convection velocity, 0.57U^. Since the space between 
neighboring microphones is 1.2hf, the phase angle 
difference between the two microphones would be, in 
the representative wavelength case, roughly 0.14 tl 
F lence, there was no need to unwrap the angles. 

Figure 12 displays a plot of the phase angle 
obtained using pairs of neighboring microphones as a 
function of the distance along the splitter plate for five 
different frequencies. By adding the phase angles from 
each of the pairs of microphones moving in the 
downstream direction, the phase shift relative to 
microphone #0 (closest to separation) was determined. 
The phase angles in the plot are calculated within the 
high-coherence frequency range. The distance along 
the splitter plate is referenced to the reattachment 
location. 



Figure 12. Streamwise development of the phase angle 
for five different frequencies (reference microphone at 
(x- x r )/x r = -0.98) 
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Except for the small region (x-x r )/x r < -0.8 in 
Figure 12, the phase plot reveals a steady, linear-like 
increase of the phase angle in the downstream direction 
for the four largest frequencies. The slopes of these 
lines differ, depending on the frequency. This indicates 
that there is a difference in the convective velocities 
corresponding to different frequencies. The convective 
velocity can be calculated by determining the slope of 
the phase plot line and using the following equation 


U 


c 


27rfx r 

AO/ 

/Ax 


(4) 


where A0/Ax is the slope deduced from the phase plot 
using a least-squares line fit. A couple of points 
regarding this calculation should be noted here: (1) the 
line fit was restricted to about (x-x r )/x r > -0.6 where the 
phase characteristics behaved linearly, and (2) because 
the horizontal axis in Figure 12 is normalized by x r , the 
calculated slope was divided by x r to convert to 
physical units. The convection velocities for the 
different frequencies are given in Table 1. Averaging 
the five convection velocities give 0.64LL, which 
differs from the convection velocity 0.57LL determined 
from the cross-correlation. This difference could be 
due to the fact that the velocities used in the average are 
not weighted by the relative dominance (energy) of the 
associated pressure fluctuations. 


Table 1. Convective velocities for various frequencies 


f(2H)/LL 

U c /U„, 

0.03 

0.63 

0.17 

0.58 

0.23 

0.66 

0.25 

0.66 

0.29 

0.68 


At f(2H)/LL = 0.17 in Figure 12, the phase 
drops below zero. This is indicative of an upstream 
convecting velocity. Fleenan and Morrison 4 is the only 
other investigation reporting negative phase angles, in 
their study of a backward-facing step. Lee and Sung^ 
mention the idea of an upstream convecting velocity but 
found no evidence of this phenomenon in their study. 
The largest negative phase angle at f(2Fl)/LL = 0.17 is 
found around microphone #4. After further 
investigation at various frequencies (ranging from 
f(2H)/LL = 0.07 to f(2H)/LL = 0.22), it was seen that 
the largest negative peak fluctuated between 
microphones #4, 5, and 6. These microphones are 


roughly around a 14 of the way to reattachment. 
Interestingly, this region was also identified earlier in 
both the RMS and the auto-correlation results, upstream 
of the changeover region. 

In the case of f(2FI)/LL = 0.03, which 
corresponds to the low-frequency peak identified in the 
power spectra, the phase angle drops below zero but has 
a maximum negative peak farther downstream of 
microphones #4, 5 and 6. The peak is located near 
microphone #11. This indicates that the upstream 
convection velocity seen in the region spanning from 
the fence to microphone #4 at frequency f(2H)AJ oc = 
0.17 extends farther downstream to microphone # 1 1 at 
frequency f(2H)/LL = 0.03. Therefore, flow structures 
at very low frequencies can be seen, as evidenced by 
the phase plot, convecting upstream starting from a 
distance x/x r = 0.53 downstream of separation. It is not 
clear if this could be associated with the flapping of the 
shear layer, which has been hypothesized to correspond 
to the low-frequency peak in the spectrum by various 
authors including Castro and Haque 6 . 

To explore the convective characteristics of 
the surface pressure around x/x r = 14 further, the cross- 
correlation function for all 28 microphones relative to 
the middle microphone #5 (which is located at about 
x/x r = 0.26) was calculated and plotted in Figure 13. 
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(X-X r )/X r 

Figure 13. Cross-correlation-coefficient contour plot 
for all 28 (reference microphone at (x-x r )/x r = -0.74) 

From the contour plot, an inclined positive 
peak is found on either side of the high correlation 
coefficient peak at microphone # 5 . These two peaks 
have opposite signed slopes with respect to each other. 
This is evidence that there are two convecting 
velocities: one upstream and one downstream as also 
deduced by Fleenan and Morrison 4 from phase 
measurements. The present study is the first to depict 
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the upstream convection from the space-time 
correlation function. Additionally, the evidence of 
downstream-traveling motion in Figure 13 shows that 
the flow structures dominant farther downstream from 
microphone #5 (those corresponding to the power 
spectrum peak at f(2H)/U„, = 0.1-0.15) are detectable as 
early as microphone #5. To check for the earliest 
manifestation of the downstream convecting motion, 
the cross-correlation map is obtained for all 28 
microphones relative to microphone #0 as shown in 
Figure 14. The axes are the same as in Figure 13. 

The cross-correlation results in Figure 14 
reveal that the downstream convecting velocity begins 
roughly around the !4x r distance, as evidenced by the 
negative-inclined contour that starts at this position. 
Upstream of 14x r , there are no negative-inclined lobes 
indicating a downstream convecting velocity. 

Therefore, it is reasoned that the flow structures seen to 
dominate the measurements downstream are first 
noticeable in the surface pressure measurements around 
the VaXj distance. 
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Figure 14. Cross-correlation-coefficient contour plot 
for all 28 (reference microphone at (x-x r )/x r = -0.98) 

The upstream convecting velocity was 
determined to be 21% of the free stream velocity. This 
is similar to the value determined by Fleenan and 
Morrison. Their value was about 20% of the free 
stream velocity. The downstream convecting velocity 
was calculated to be 47% of U„. This value is lower 
than 57% U„, calculated earlier from the cross- 
correlation map with the reference microphone located 
closest to reattachnrent (Figure 10). Fleenan and 
Morrison 4 also noticed this trend in their impermeable 
backward-facing step. Upstream near the fence, they 
observed a 0.5U*. convecting velocity that increased to 


0.6U„, close to reattachment and continued to rise 
farther downstream. 

Concluding Remarks 

The space-time characteristics of the surface- 
pressure within the separating/reattaching flow region 
of a splitter-plate-with-fence configuration were 
studied. For this purpose, a comprehensive database 
was compiled using an 80-microphone array embedded 
in the wall of the splitter plate. Only the time-averaged 
space-time statistics of the surface pressure 
measurements have been presented in this paper. In 
general, the results from this analysis compared well 
with available literature in related, but not exactly 
similar flow geometries. 

The streamwise distribution of RMS pressure 
fluctuations exhibited a rapid rise in magnitude around 
‘/4X r , with the peak value occurring in the vicinity of 
reattachment. The region extending up to l Ax r was also 
identified in the auto-correlation function analysis, 
which showed decreasing time scales with downstream 
distance. The R p - p - contour plot revealed that within the 
'Ax r region the surface-pressure signature was 
dominated by large time scales. Farther downstream, 
near reattachment, smaller time scales were prominent 
in the wall-pressure measurements. Transition between 
the two different time scales occurred in a region 
extending from 0.25x r to 0.7x r . 

The peak energy in the power spectra for 
microphones close to the fence is concentrated at very 
low-frequencies (f(2H)/U 00 = 0.02 - 0.03), which has 
been attributed by many researchers to the ‘flapping’ of 
the shear layer. Farther downstream near the mean 
reattachment location, the concentration of peak energy 
in the power spectra is seen at a higher frequency 
(f(2H)/U«, = 0.1-0.15), relating to the highly turbulent 
structures within the shear layer. The findings from the 
auto-correlation function and power spectra analyses 
were consistent. 

From the cross-correlation function analysis, 
with respect to a reference microphone located closest 
to reattachment, the wall-pressure signature of the 
downstream convective motion of shear-layer 
structures, described by Cherry et al} amongst others, 
was found to travel at a convection velocity of 0.57U„. 
Phase angle information, obtained from the cross- 
spectrum of the signal between neighboring pairs of 
microphones and constrained to frequencies up to 
f(2H) AJ^ = 0.3 (range of high signal coherence), 
revealed an upstream convecting velocity of 0.2 1U M at 
low frequencies within the 'Ax r region. A cross- 
correlation contour plot, based on a reference 
microphone located near l Ax r , showed the two opposing 
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convection velocities and provided evidence that the 
earliest detection of the downstream-traveling turbulent 
structures seen near mean reattachment was around 

'/4X r . 

Overall, two distinctive regions emerged from 
the spatio-temporal analysis of the surface pressure 
measurements in both the time and frequency domains. 
Upstream, from the fence to 'A x r , the surface-pressure 
signature was dominated by large time scale 
disturbances and an upstream convecting velocity of 
0.21U^. Beyond the 1 Ax r , turbulent structures with 
small time scales and a downstream convection velocity 
of 0.57U^ generated most of the pressure fluctuations. 
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